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We study the decay of entanglement of quantum dot electron-spin qubits under hyperfine interac- 
tion mediated decoherence. We show that the system displays two qualitatively different behaviors 
in the low and high magnetic field limits. At high magnetic fields, an initial Bell state will undergo 
exponential entanglement decay, where the amount of entanglement left in the system is propor- 
tional to the non-zero coherence. Contrarily, at low magnetic fields sudden-death-type behavior is 
observed and disentanglement is faster than the decay of coherence. Furthermore, we have found 
that the time of complete disentanglement displays a counter-intuitive, non-monotonic dependence 
on the magnetic field. We point out a directly measurable entanglement witness which is in this 
case sensitive to the magnetic field, and can serve as a threshold sensor for the magnetic field. 
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Systems of electron spins confined in quantum dots 
(QDs) have received much theoretical (see Refs [I-Q for 
review) and experimental [J-ul interest since the initial 
proposal for spin-based quantum computing [8| . It is now 
known that the primary source of decoherence in such 
systems is the hyperfine interaction between the electron 
spin and the spins of the nuclei of the QD atoms. This 
interaction leads to pure dephasing at moderately high 
magnetic fields, but at lower magnetic fields a more in- 
volved decoherence process is seen which leads to a re- 
distribution of the electron spin occupations. 

At low magnetic fields the calculation of the evolu- 
tion of the electron spin becomes difficult, because of 
the "flip flop" terms in the hyperfine-interaction Hamil- 
tonian. Now these terms cannot be treated perturba- 
tively (the interaction can be regarded as a small pertur- 
bation with respect to the electron Zccman splitting only 
at moderately high magnetic fields [2,^]). Furthermore, 
the dependence of the hyperfine coupling constants on 
the location of a given nucleus with respect to the elec- 
tron wave function substantially complicates the prob- 
lem. This combined with the discrete structure of the 
nuclear spin environment [9[ makes the problem strongly 
involved numerically. 

Due to such complications, the study of spin entangle- 
ment [l3| decay in a two electron two QD system, which 
is essential for any postulated qubit for quantum compu- 
tation schemes, has been limited to a number of complex, 
yet solvable scenarios [ll| - [l3j . The complexity accounts 
for the nontrivial behavior of the reported evolution of 
entanglement. Hence, in Refs [11| and [Tj] the uniform 
coupling ("box") model is extended to account for the 
exchange interaction between electron spins for a small 
number of nuclei in the common nuclear bath limit (with 
low bath polarization) and separate nuclear baths limit 
(with high bath polarization and large exchange inter- 
action), respectively. Ref. [ij] utilizes the "box" model 



with a simplified thermal spin bath state to introduce a 
scheme for multipartite entanglement generation medi- 
ated by the interaction with a nuclear bath. An excep- 
tion is Ref. [Ijl , where the evolution of entanglement of 
non-interacting spin qubits is studied, but the decoher- 
ence model considered is phenomenological and leads to 
a different type of decoherence than is reported in the 
literature [2|, Q ■ The importance of Ref. [ij] lies in its 
attempt to quantify multipartite entanglement under a 
feasibly realistic evolution of the qubit states. 

We study the evolution of entanglement of two non- 
interacting electron spin qubits confined in two well sep- 
arated lateral GaAs QDs. The qubits interact with sepa- 
rate nuclear spin reservoirs via the hyperfine coupling, 
which are taken in the high-temperature thermalized 
state to which the baths relax quickly at experimentally 
accessible temperatures [^]. In this case, we can use the 
"box" model for the whole range of magnetic field val- 
ues |16l . |l7[. This is because the hyperfine- interaction- 
induced entanglement decay takes place on time scales 
much shorter than the "box" model limit of applicabil- 
ity, t < N/A, where N denotes the number of nuclei, and 
A = ^j, Ak is the sum of coupling constants between the 
electron and the nuclei. 

We show that the nature of the decay of entangle- 
ment changes substantially when the transition to the 
high magnetic field limit is made. To this end, we study 
entanglement decay of an initial Bell state, for which sud- 
den death of entanglement [i3,[i3 (complete disentangle- 
ment while the loss of coherence is still only partial) is not 
possible under pure dephasing processes [20| . In the high 
magnetic field regime, where decoherence is restricted to 
pure dephasing, entanglement is proportional to the co- 
herence as expected, and decays following the appropri- 
ate exponential function. Contrarily, at low magnetic 
fields the evolution of the electron spin state involves a 
redistribution of the occupation of the spin-up and spin- 



down states. Although the decay of single qubit coher- 
ence is qualitatively similar in both regimes when the 
spin-up and spin-down occupations are equal, the result- 
ing low-magnetic-iield Bell state disentanglement is very 
different from the high-magnetic-field disentanglement, 
and leads to the appearance of sudden death. Hence, the 
same system realizes qualitatively different disentangle- 
ment scenarios in different magnetic field regimes for the 
same initial entangled state. 

The system under study can be described by a sepa- 
rable Hamiltonian, H = Hi (E)l2 + Ii ^ H2, where the 
individual QD subsystems are described by Hamiltoni- 
ans of the form (the magnetic field is applied in the z 
direction) , 

H, = -gt,BS!B + Y,Ak,.S^il, (!) 

k 
k 

with the index i = 1, 2 distinguishing the two dots. The 
first term in the Hamiltonian ([T]) is the electron Zeeman 
splitting, where g is the effective electron g-factor, /is is 
the Bohr magneton, S^ is the component of the electron 
spin parallel to the magnetic field, and B denotes the 
applied magnetic field. The last two terms describe the 
hyperfine interaction between the spin of an electron and 
the spins of the surrounding QD nuclei. The diagonal 
(second) term is also known as the Overhauser term and 
leads to pure dephasing, while the last term, known as 
the "flip-flop" term, is responsible for both dephasing and 
leveling out of the electron spin occupations. Here, I^.i 
are spin operators of the individual nuclei (discriminated 
by the index k) in dot i. /| j is the component parallel to 

the applied magnetic fleld, while I^. .^ = I^j^±ill,^ are the 
nuclear spin raising and lowering operators. Analogously, 
S'j = Sf ± iSf arc the raising and lowering operators for 
the electron spin. The coupling constants of the hyperfine 
interaction depend on the species of the nuclei and on its 
location with respect to the electron wave function. 
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where A^ ^ = |/xo7e7fc,i are the coupling constants of 
a given nuclear species found at site k of dot i, with 
liQ denoting the vacuum magnetic permeability, 7e and 
jkji being the electron and nuclear gyro-magnetic ratios, 
respectively, while vq is the unit cell volume of the QD 
crystal, ^^ (r) is the wave function of the electron located 
in dot i, and Tk.i is the position of the fc-th nucleus in 
dot i. 

We have omitted the nuclear Zeeman term and the 
dipolar interaction between nuclei in the Hamiltonian ([T|) . 
The first was omitted because nuclear Zeeman energies 
of gallium and arsenic are very small, and the resulting 
energy splittings are of the order of tens of neV (corre- 
sponding to less than a mK) for each Tesla of magnetic 
field applied to the system. The nearest neighbor dipolar 



coupling constants between nuclei are even smaller, and 
are of the order of 0.1 neV. Hence, at typical experimen- 
tal temperatures both nuclear terms in the Hamiltonian 
are much smaller than kgT [3l. Il6l|. 

For the same reason the nuclear baths can be de- 
scribed by infinite-temperature, fully mixed density ma- 
trices [3|, |9| unless the state of the nuclear environment 
is especially experimentally prepared. While a polarized 
environment strongly changes the resulting dynamics and 
leads to an increase of the electron spin coherence time 
jlSl l2l| . the preparation of such an environment is de- 
manding experimentally J22ll23l | and the currently attain- 
able levels of polarization are under 70% [23| ■ The study 
of polarized environments is beyond the scope of this pa- 
per, where we wish to describe spin disentanglement in 
the simplest and most common scenario. Hence, we limit 
the study to initial states, for which the density matrix of 
the two-qubit subsystem and the two nuclear reservoirs 
is in the product state q{0) = pdqd{0) ^ -Ri(O) (8) i?2(0), 
where Pdqd{0) is the initial state of the two confined 
electron spins, while the nuclear baths Ri are initially 
fully mixed. 

In the following, we use parameters corresponding to 
two identical lateral GaAs QDs, but the results are qual- 
itatively valid for any dot type, as long as they can 
be treated as non-interacting. Electron wave functions, 
which are necessary to find the coupling constants of the 
hyperfine interaction, Eq. ^ , are modeled by anisotropic 
Gaussians with the extension ^j^ = 20 nm in the xy plane 
and /z = 2 nm along z direction, which is the direction of 
the applied magnetic field. The number of crystal unit 
cells considered within each dot is iVi = A'2 ~ 1.5 x 10^. 

All isotopes naturally found in GaAs carry spin / = 
3/2 and the nuclear-species-dependent coefficients A°j. 
are equal to Aq^gb = 36 /xeV, Aq^l''^ = 46 peY, and 
AasTs = 43 /xeV [H, [13 for the different types of nu- 
clei. The relative abundances of the gallium isotopes are 
60.4% for Ga'^^ and 39,6% for Ga'^^ , this together with 
the fact that there is one gallium and one arsenic atom in 
the GaAs unit cell gives the average hyperfine coupling 
constant A = 83 fieV. The g-factor is equal to g = —0.44 
[23, hence, the Zeeman electron spin splitting is equal to 
25.5 ^eV per Tesla of magnetic field. 

The parameters are used to find single QD evolutions 
in the high-magnetic-field limit, gpsB > A, for which the 
"fiip-flop" term may be completely neglected. The condi- 
tion is fulfilled for magnetic fields greater than about 3.25 
T. The Hamiltonian is then diagonal and it is possible to 
find the QD state evolution for a realistic distribution of 
coupling constants Ak given by Eq. ^ while taking into 
account the large number of nuclei N ^ 10^. The result- 
ing dynamics is limited to pure dephasing which is fur- 
ther independent of the magnetic field (and local unitary 
oscillations that are irrelevant for the study of entangle- 
ment) for the initial high-temperature density matrix of 
the nuclear baths. As predicted [la], the decay of a single 
spin is proportional to exp(— i^/T2*^), with a characteris- 
tic time constant T,* — ^ / j/j.^ 



N/A, where the same 



spin / for all nuclei is assumed. vN/A « 10 ns accord- 
ing to the parameters introduced in the last paragraph 
and the Tj* = 12.36 ns extracted from the calculation 
corroborates this. 

To quantify the single dot evolution at lower mag- 
netic fields, we use the "box" model which is valid on 
short time scales when the high-temperature nuclear bath 
density matrix is used and at high magnetic fields con- 
verges with the approach introduced in the last para- 
graph. The upper limit of short-time-scale behavior is 
approximated by N/A [23|, the value of which is 1.2-10^ 
ns for the parameters used and exceeds the disentangle- 
ment times by three orders of magnitude. In the "box" 
model, the hyperfine coupling terms are assumed con- 
stant Ak ^ a = A/N, which allows for exact diagonal- 
ization of the Hamiltonian ([ij as outlined in the Sup- 
plementary Materials. Furthermore, "box" model evo- 
lutions involving a large number of nuclei can be suc- 
cessfully simulated with reasonably small numbers of nu- 
clei, since few-body coherent effects disappear already in 
the case of 10 spins-3/2 and for 50 spins large-number- 
of-nuclei evolutions are reproduced (see Supplementary 
Materials for details). 

The type of evolution manifested by the single QD sys- 
tem depends strongly on the value of the magnetic field. 
At very low magnetic fields, QD occupations are partially 
leveled out due to the interaction with the environment. 
In this regime, phase decoherence closely resembles the 
decay of the occupation difference. The effect of the en- 
vironment on the occupations is diminished with growing 
magnetic field, while coherence damping remains strong, 
although it starts to resemble exponential decay. When 
the limit of high magnetic fields is reached, the interac- 
tion with nuclear spins cannot disturb the QD occupa- 
tions, and the pure dephasing process follows a Gaussian 
decay proportional to exp(— t^/T2*^). 

The study of the evolution of entanglement of the two 
QD system requires a two-qubit entanglement measure 
which can be calculated from the system state. One such 
measure, for which an explicit formula is available, is the 
concurrence [2^, |30l , which is closely related to the en- 
tanglement of formation, defined as the ensemble average 
of the von Neumann entropy minimized over all ensemble 
preparations of the state [Sll, [32| . The concurrence for bi- 
partite entanglement is given by C{pdqd) = max{0, Ai — 
A2 — A3 — A4}, where A^ are the square roots of the eigen- 
values of the matrix PDQD{(^y®<yy)p*DQD{^y®'^y)- Here, 
Pdqd is the two qubit density matrix, P^jqd i^ i^^ com- 
plex conjugate, and ay is the appropriate Pauli matrix. 

We will quantify the entanglement evolution of initial 
maximally entangled Bell states, |*^) = l/%/2(|l) ± |2)) 
and 1$*) = 1/a/2(|0)±|3)), where the states in the single 
QD basis are equal to |0) = |tt), |1) = |ti), |2) = lit), 
and |3) = I H), under the influence of the interaction 
with separable nuclear reservoirs in the high tempera- 
ture regime. The evolution of the coherences for these 
initial states is limited to the single off-diagonal element 
of the density matrix which is initially non-zero, while 
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FIG. 1: Time evolution of entanglement for different magnetic 
field values, B = - solid red line, B = 0.0525 T - long 
dashed green line, B — 0.08 T - dashed blue line, B — 0.095 
T - dotted pink line, and B = 1 T - dashed/dotted blue line 
(high magnetic field limit). 



the other coherences remain zero at all times. Contrar- 
ily, all four occupations are influenced (except for the 
high-magnetic-field limit where the decoherence is a pure 
dephasing process) by the interaction. Hence, the double 
QD density matrix is simplified and the concurrence is 
always of the form, 



C{pdqd) = 2max{0, |p.y | - ^/pkkPii}, 



(3) 



where i,j are equal to 1, 2 or 0, 3 depending on the ini- 
tial state, and k y^ I, k ^ i, k ^ j, I y^ i, I ^ j. It 
is evident from Eq. ^ that sudden death of entangle- 
ment will occur when \pij\ < ^PkkPii, so it is expected in 
the low magnetic field regime when the QD occupations 
are disturbed, while it will not occur for high magnetic 
field pure dephasing. Furthermore, because the qubits 
interact with separate environments at high temperature 
thermal equilibrium, the evolution of entanglement is the 
same for all four initial Bell states. 

Fig.[T]shows entanglement decay for different magnetic 
field values. The zero magnetic field curve (red solid) 
limits all higher magnetic field curves from below and 
ends in sudden death. The high magnetic field curve 
(dashed/dotted blue line), provides the upper limit for 
the concurrence at a given time and undergoes expo- 
nential decay. In between, the curves corresponding to 
small magnetic fields display a more complex entangle- 
ment evolution. According to ([3]), the visible oscillations 
arc due to the interplay of the dephasing process and the 
shifts in the occupations. The number of oscillations in- 
creases with the increase of the magnetic field, while they 
become less pronounced in the process. This is because 
high magnetic field prohibits changes in occupation lev- 
els. Due to the oscillations, the sudden death times in 
this regime are not a monotonous function of the mag- 
netic field. 
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FIG. 2: Entanglement sudden death time tsD as a function 
of the magnetic field. The vertical lines mark the magnetic 
field values corresponding to the long dashed green, dashed 
blue and dotted pink lines in Fig. [T] 



This is quantified in Fig. [51 where the dependence of 
the sudden death time is plotted as a function of the 
magnetic field (the magnetic field values corresponding 
to the long dashed green, dashed blue and dotted pink 
lines in Fig. [T]arc marked by the vertical lines). At low 
magnetic fields, an strong oscillatory behavior is evident, 
starting from around 10 mT. The oscillations disappear 
at around 95 mT, long before the high magnetic field 
limit is reached, and the times of sudden death increase 
quickly and almost monotonously with a further increase 
of the magnetic field, until the pure dephasing (no sudden 
death) limit is reached. 

For any fixed time the considered dynamics is a ten- 
sor product of two quantum channels (completely posi- 
tive trace preserving maps) which are both bi-stochastic 
- they preserve the maximally mixed state. The proof 
of this statement is straightforward, namely: Consider a 
single qubit subsystem in a maximally mixed state cou- 
pled by an arbitrary unitary interaction to the maximally 
mixed environment state. Since the product of identities 
is the identity operator on the composed (qubit -I- en- 
vironment) system and is invariant under any unitary 
operation, the final state of the two subsystems remains 
the same product of maximally mixed states. In par- 
ticular the maximally mixed state of the qubit subsys- 
tem is preserved. The above argument applies to each 
of the QDs separately. It follows that any two-qubit 
state with maximally mixed subsystems will retain the 
property during the evolution. A system with maxi- 
mally mixed subsystems has to have Bell diagonal qubit 
states [33|, and starting from the Bell diagonal state 
guarantees Bell diagonality for the whole evolution. It 



is known that when a Bell diagonal state becomes sep- 
arable in the course of its evolution, the entanglement 
fidelity F becomes 1/2. We start from a single Bell 
state (a trivial Bell diagonal state) and the maximal 
eigenvalue corresponds to the projection onto that state. 
Monitoring the difference of that eigenvalue and 1/2, i.e. 
W{t) = F - \ ^ (*o|p(i)l*o) - 1/2 (which is a spe- 
cific entanglement witness), we may identify the moment 
of entanglement sudden death exactly. This means that 
the zero point t* of W{t) is just the sudden death time, 
t* = tsD, and as such has the same dependence on the 
magnetic field as shown in the Fig. [5]). Quite remarkably, 
W as an entanglement witness is directly measurable. By 
measuring this quantity we can get the exact estimate of 
the magnetic field whenever it corresponds to the initial 
monotonic regime of the function. In the regions close to 
the steep parts of the function, the above value is quite 
sensitive to the field B and can be considered as a thresh- 
old sensor of the magnetic field. 

We have studied decay of QD spin-qubit Bell state 
entanglement under decoherence processes mediated by 
the hyperfine interaction between confined electron spins 
and the spins of the QD nuclei in the limit of separate 
nuclear reservoirs and non-interacting QDs. We have 
shown that varying the magnetic field leads to a transi- 
tion between substantially different entanglement decay 
processes; sudden death of entanglement is manifested 
at low magnetic fields, but at sufficiently high magnetic 
fields exponential entanglement decay is observed. Fur- 
thermore, at low magnetic fields, the evolution of en- 
tanglement displays counter-intuitive oscillatory behav- 
ior which results in a non-monotonic dependence of the 
sudden death time as a function of the magnetic field. 
The characteristic behavior is an outcome of the interplay 
of the decay of the system coherence with decoherence in- 
duced redistribution of the double QD spin occupations. 
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In the "box model" , the hypcrfinc couphng terms are 
assumed constant Ak — a ~ ^/N, which allows for the 
exact diagonalization of the single dot Hamiltonian (Eq. 1 
in the main article) as outlined below. 

Firstly, it is now possible to rewrite the single dot 
Hamiltonian in terms of the components of the total nu- 
clear spin operator Ki = ^^, /^ ,j. Furthermore, it is con- 
venient to use the eigenstates of the total nuclear spin and 
its z-component as the nuclear environment basis states, 
WK, to)}, described by the spin quantum numbers K and 
TO., TO. = —K,—K-^\,...,K, and fulfilling the relations 
k'l\Km) = h^K{K + l)\Km) and Kf\Km) = hm\Km). 
Note that the Hamiltonian (Eq. 1 in the main article) 
acting on any state \a\K,m), where a =tj4- denotes the 
electron spin, cannot change the nuclear quantum num- 
ber K and conserves the z-componcnt of the total spin of 
the combined electron and nuclear spin system. Hence, 
the Hamiltonian can be represented in easily diagonaliz- 
able 2x2 block form, where each block links the states 
I \]K,m) and | J,; iiT, to, + 1), which form a closed sub- 
space for every K and to e [—K,K — 1]. The form of 
these 2x2 blocks is given by 



Em MK,m 

MK,m —Ejn+l 



where the energies are given by £",„ = h/2{Q, + 
am), with the Zeeman frequency f2 — —gfj,BB/h, 
and the transitions are governed by Mx^m ~ 

The eigenvectors are 



na/2.jK{K + 1) - to(to + 1). 
then of the form 



\+;K,m) = cos6'_ft:,m| 1;;K,m) +s\\i9K,'m\ i;K,Tn + 1), 
\-;K,m) = -sm9K,m\ \]K,m) +cos6'K,m| i; A',to + 1), 



with 



M 



sm ok.i 



Ka 



(Etm+Em+ir+MlJ 

and with corresponding eigenvalues given by 



(1) 



E 



ha/2±J{Em + E,rr+iY+^M\ 



K,m 



any combined state of the electron and nuclear spins can 
be found. 

The high-temperature nuclear spin density matrix, 
rewritten in the {\K,m)} basis takes the diagonal form 



R{Q)=Y,PK,m\K,m){K, 



(3) 



Ka 



where PK,m are coefficients describing the multiplicity 
of the occupation of each state and satisfy the relation 
'^Kjn^K-^n ^ 1- ^^^ nuclei with spin s they are given 
by|f 



^^,^--^(-1)' 



N\ ({s + l)N - (2s + l)i - A" - 2 

N -2 



(4) 
where i G [0, N\ is an integer. For spin i systems the 
formula simplifies to 



Pif,, 



N\{2K+l) 



{\N - K)\aN + K + l)\' 



(5) 



The evolution of the QD density matrix can be found 
by tracing out the environmental degrees of freedom and 
assuming a product initial state p(0) = /oqd(O) ® i?(0), 
with the initial QD density matrix denoted by pq_d(0) 
and the nuclear spin density matrix i?(0) given by eq. 
^. It is described by 

PQD{t) = |t>(tl(m(o)^W + Pu(o)(i-^'(t))) 

+U)ai(pn(o)(i-^(t)) + Pu(o)^'W) 
+\^){i\{Pn{mit)) + H.c.. (6) 

Here, X{t) = Y.K,nPK,m\XK,m{t)\\ 

X'{t) = Y.K.mPK,m\XK^m-i{t)\'', and Y{t) = 

J2K,rnPK:^nXK,mit)XK,m^l{t), with 

XK,rait) = cos^ e^.^e-^^^'-* + siu^ 0/^,™e * *^'."*. 
The phase coeflficients are equal to 



K,m 



(2) 



')K,m — 



hJn^ + na{2m + 1) + ^ + a^K{K + 1) 



Furthermore, the | f; K, K) and | |; K, —K) states are 
eigenvectors of the Hamiltonian. Hence, the evolution of 



and the mixing angles 9K,m are given by eq. ([T]). Note 
that for a large number of nuclei N the time-dependent 
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FIG. 1: "Box model" evolution of the difference of single QD 
occupations Ap(t)/Ap(0) (see text) for B = mT (a), 1.5 mT 
(b), 15 mT (c), and 150 mT (d). Different curves correspond 
to different numbers of spin-3/2 environment nuclei, N = 2 
(red, solid line), 10 (green, dashed), and 50 (blue, dotted). 
Points denote the evolution fitted with 300 nuclei with spin 
1/2. 



FIG. 2: "Box model" evolution of the amplitude of single QD 
coherence Ptiit) / Pfi{0) (see text) for 5 = mT (a), 1.5 mT 
(b), 15 mT (c), and 150 mT (d). Different curves correspond 
to different numbers of spin-3/2 environment nuclei, N = 2 
(red, solid line), 10 (green, dashed), and 50 (blue, dotted). 
Points denote the evolution fitted with 300 nuclei with spin 
1/2. 



functions X{t) and X'{t) coincide, and the evolution of 
the diagonal terms of the density matrix is described by 



Ap(t) = Ap(0)(2X(t)-l), 



(7) 



with Ap{t) = ptt(0 - Pliit)- 

Even though the "box model" is exactly solvable as 
reproduced above, finding the actual QD evolution when 
the environment is in the high-temperature equilibrium 
state becomes numerically challenging very quickly with 
the growing number of nuclei N, due to the involved 
summation over K and m, and is practically impossi- 
ble for a realistically large values of N. The difficulty 
of the task also grows rapidly with higher spins of the 
nuclear species taken into account. Figs. [T] and [2] serve 
to demonstrate that "box model" evolutions involving 
a large number of nuclei can be successfully simulated 
with reasonably small numbers of nuclei, since the few- 
body coherent effects disappear already in the case of 
10 spins-3/2 and for 50 spins the large number of nu- 
clei evolutions are reproduced. The necessary condition 
to achieve convergence is that A/^/N remains constant. 
This requirement stems from the semi-classical approx- 
imation [2| result giving the characteristic decay time 

N /A, where I is the nuclear spin. 



T. 



6 



^(^+1) 



Fig. [T] contains plots of the function 1 — 2X{t) which 
determines the evolution of the QD occupations for large 
N (see eq. (O) for four different magnetic field values 
as a function of time. We have found that the function 
X'{t) converges for similar values of N as the function 
X{t), hence no additional plots are necessary. Analo- 
gously, the plots of the evolution of the amplitude of the 
off-diagonal terms, \Y{t)\ = |pti(0/Pti(0)l ^^^"^ shown in 



Fig. [2] for the same values of the magnetic field. In each 
plot, there are three curves corresponding to spin-3/2 
nuclei, with A^ = 2, 10 and 50. We have also found the 
evolutions for N = 500 and N = 2000, but the resulting 
curves are indistinguishable from the 50 nuclei curves on 
relevant time scales and, hence, have not been included 
in the plots. As can be seen, the effects resulting from 
a limited number of spins are strong only for very small 
N. They manifest themselves as additional oscillations 
of the 2 nuclei occupation curves and echo-like charac- 
teristics resulting from the alignment of the few nuclear 
spins at certain time intervals. 

The type of evolution manifested by the single QD sys- 
tem depends strongly on the value of the magnetic field. 
At very low magnetic fields, the occupations of the QD 
are redistributed due to the interaction with the environ- 
ment, see Fig. [T] (a) and (b). In this regime, the phase 
decoherence closely resembles the decay of the occupation 
difference. Fig. [2] (a), (b). The effect of the environment 
on the occupations is diminished with growing magnetic 
field. Fig. [T] (c), (d), while coherence damping remains 
strong, although it starts to resemble exponential decay. 
Fig. [2] (c), (d). When the limit of high magnetic fields 
is reached, the interaction with nuclear spins cannot dis- 
turb the QD occupations, and the pure dephasing process 
follows a Gaussian decay proportional to exp(— t^/Tj*^). 

Analysis of spin- 1/2 environments further simplifies 
the generation of QD system evolutions for a given num- 
ber of nuclei. To quantify the applicability of such an 
approximation, points have been added to the plots of 
Figs. [T] and m which denote the evolutions found by mod- 
eling the spin 3/2 environment with a spin 1/2 nuclei for 
N = 300. The large number of nuclei is necessary to 



achieve convergence for spins 1/2. The fitting required a 
scahng of the constants a by \ -r^^-Tr^^-rn = VS, with 
I3.I') — i, I1/2 — 55 deduced from semi-classical approx- 



A/2(-fl/2 + l) 

3/2 - 2' ^1/2 

imation. The resulting evolutions are qualitatively and 
quantitatively reproduced very well and the transitions 
between different types of decoherence with growing mag- 



netic field are the same as in the case of spin-3/2 envi- 
ronment. Unfortunately, the large number of environ- 
ment atoms required diminishes any computational ad- 
vantages, which would be gained by using nuclei with 
smaller spins. 



[1] V. V. Miliailov, Journal of Physics A: Mathematical and 

General 10, 147 (1977). 
[2] I. A. Merkulov, A. L. Efros, and M. Rosen, Phys. Rev. B 



65, 205309 (2002). 



